Abstract. In a previous study (Maurin et al., 2001), wc explored the set of parameters 
describing diffusive propagation of cosmic rays (galactic convection, reacceleration, halo 
thickness, spectral index and normalization of the diffusion coefficient), and we identified 
those giving a good fit to the measured B/C ratio. This study is now extended to take 
into account a sixth free parameter, namely the spectral index of sources. We use an 
updated version of our code where the reacceleration term comes from standard minimal 
£>\j ' reacceleration models. The goal of this paper is to present a general view of the evolution 
, of the goodness of fit to B/C data with the propagation parameters. In particular, we 
, find that, unlike the well accepted picture, and in accordance with our previous study, 
i^P a Kolmogorov-likc power spectrum for diffusion is strongly disfavored. Rather, the x 2 
analysis points towards S > 0.7 along with source spectra index < 2.0. Two distinct 
, energy dependences are used for the source spectra: the usual power-law in rigidity and 
a law modified at low energy, the second choice being only slightly preferred. We also 

m ■ 

J> • show that the results are not much affected by a different choice for the diffusion scheme. 
qq ' Finally, we compare our findings to recent works, using other propagation models. This 
| study will be further refined in a companion paper, focusing on the fluxes of cosmic ray 
, nuclei. 
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1. Introduction 

Cosmic rays detected on Earth with kinetic energies per nucleon from 100 MeV/nuc to 
100 GeV/nuc were most probably produced by the acceleration of a low energy galactic 
population of nuclei, followed by diffusion in the turbulent magnetic field. The accelera- 
tion process and the diffusion process have a magnetic origin, so that they should depend 
on rigidity. The rigidity dependence of the diffusion coefficient is given by quasi-linear 
theory as 



where the parameters Kq and S should ideally be given by the small-scale structure of the 
magnetic field responsible for the diffusion. As this structure is not well observed, some 
theoretical assumptions must be made in order to predict 6. As regards the spectrum 
just after acceleration, the situation is far from clear, as it depends on the details of the 
acceleration process. Several models give a power-law distribution (e.g. Bcrczhko et al., 
1994, Gieseler et al., 2000) 

^ oc TZ- a (2) 
dp 

with a definite value for a which depends on the model. 

Most analyses of cosmic ray nuclei data assume given power-laws for the diffusion and 
acceleration energy dependence , so that the results partially reflect certain theoretical 
a priori. In this work, we try to avoid this bias by determining the quantities a and S 
directly from the data, in particular B/C, for reasons exposed below. 

The paper is organized as follows. We first recall the main features of our diffusion 
model. As a few modifications have been made since previous works, §^ is devoted to 
their description and justification. Then, the analysis method is described in §0 and the 
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results are shown and discussed in a comparison is eventually made with other similar 
works in 



2. Description of the model 

This paper and its companion (Donato, Maurin & Taillet, in preparation) use the same 



description of cosmic ray propagation as our previous analyses ( Maurin et al. (2001) 



Donato ct al. (2001)| , ponato ct al. (2002)| , [Barrau et al. (2002)| , [Maurin ct al. (2002)| ). 
Particles are accelerated in a thin galactic disk, from which they diffuse in a larger vol- 
ume. When they cross the disk, they may interact with interstellar matter, which leads to 
nuclear reactions (spallations) - changing their elemental and isotopic composition - and 
to energy losses. Interaction with Alfven waves in the disk also leads to diffusive reaccel- 
cration. The reader is referred to Maurin et al. (2001) - hereafter Paper I - for all details, 
i.e. geometry and solutions of our two zone/three-dimensional diffusive model, nuclear 
parameters (nuclear grid and cross sections), energy losses terms (adiabatic, ionization 
and Coulomb losses), solar modulation scheme (force-field), as well as general description 
of the procedure involved in our fits to data (selection of a set of parameters, x 2 test 
comparison to data). In particular, though some inputs are modified (see next section), 
the final equation describing cosmic ray equilibrium is formally equivalent to that of Pa- 
per I (sec Eq. A13): it is a second order differential equation in energy solved with the 



Crank-Nicholson approach (see Ponato et al. (2001)| , Appendix B - hereafter Paper II) . 



Finally, a schematic view of our diffusion model is presented in Barrau et al. (2002) and 
Fig. ^ (see next section) summarizes the algorithm of our propagation code. 

Some aspects of this model are formally unrealistic. First, the distribution of inter- 
stellar matter has a very simple structure: it does not take into account a possible z 
distribution inside the disk (thin disc approximation is used instead), nor radial and an- 
gular dependence in the galactic plane. The orthoradial 9 dependence would even be more 
important from an accurate description of the magnetic fields and the ensuing diffusion, 
as flux tubes are likely to be present along the spiral arms. However, this is not crucial as 
we are interested in effective quantities (diffusion coefficient and interstellar density) but 
not in giving them a "microscopic" explanation. This is why we chose to use a universal 
form of the diffusion coefficient, with the same value in the whole Galaxy. Finally, it is 
known that a fully realistic model has to take into account interactions between cosmic 
ray pressure, gas and magnetic pressure, i.e. magnctohydrodynamics. 

The semi-analytical diffusion approach should be thought of as an intermediate step 
between leaky box approaches and magnctohydrodynamics simulations and is actually 
justified by these two very approaches: the first showed that the local abundances of 
charged nuclei can be roughly described by two phcnomcnological coefficients - the escape 
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length and the interstellar gas density in the box. The second hints at the fact that the 
propagation models such as the one used here arc well suited for the description of cosmic 
ray physics. 

However, it is difficult to conclude whether these parameters arc valid for other kinds 
of cosmic rays (e + , e~, nuclei induced 7-ray production) and whether they are either 
meaningful but valid only locally on a few kpc scale (i.e. not in the whole Galaxy - see 
as an illustration Brcitschwerdt, Dogiel & Volk (2002)), or meaningless but phenomcno- 
logically valid as an average description of more subtle phenomena (see as an example 



the discussion of the Alfvenic speed in Sec. 6.3.4). 



3. New settings 

Only a few ingredients differ from our previous analysis (Paper I) . The reason for these few 
changes is twofold: first, we attempt to use a better motivated form of the reacceleration 
term; second, as the real value of the exponent in the source power-law cannot be firmly 
established from acceleration models - the latter being seemingly different from what is 
naively deduced from direct spectra measurements -, it becomes a free parameter in the 
present analysis. 



3.1. Transport of cosmic rays 

The starting point of all cosmic ray data analysis is the transport equation. As emphasized 
in Berezinskii et al. (1990), a diffusion-like equation was first obtained phenomenologi- 
cally. Afterwards, the kinetic theory approach provided grounds for a consistent deriva- 
tion. This transport equation reads: 

I - V(KV, - V c f) - ^|< P V> - ^/K^f + f (3, 
In this equation, / = f(t,r,p) is the phase space distribution, K is the spatial diffusion 
coefficient, K pp is the momentum diffusion coefficient; both are related to the diffusive 
nature of the process. Finally V c is the velocity describing the convective transport of cos- 
mic rays away from the galactic plane. Actually, the full equation of cosmic ray transport 
includes other terms, such as catastrophic and spontaneous losses, secondary spallative 
contributions and continuous energy losses (coulombian and ionization losses). These 
were taken into account as described in detail in Paper I, to which the reader is referred 
for a complete description and references. They will not be further discussed here. 

This equation can be rewritten using the cosmic ray differential density dn/dE = 
N(E). As the momentum distribution function is normalized to the total cosmic ray 
number density (n = 4ir J dpp 2 f), we have N(E) = (4ir/{3)p 2 f to finally obtain 
™ - V [KVN { E) - V C N( E)] - ^QO-ipW) = (4) 
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Q(E) ; 



with 

QMccff. ^ (5) 

In this paper, K pp will be taken from the quasi- linear theory (see below). 

From a theoretical point of view, the most natural choice for the energy dependence 
of the source term seems to be a power-law in rigidity (or momentum) for dQ/dp. This 
translates into Q(E) oc R~ a / ' (3 in our set of equations (see Eqs. (|]) and (|J) above). 
Several different forms were used in the past because of the lack of strong evidence from 
observed spectra (see for example Engelmann et al. (1985), Engclmann ct al. (1990)). In 
particular, our previous analysis allowed only a rigidity dependence Q{E) oc R~ a (for 
the special case 7 = S + a as 2.8). These two forms differ only at low energy and we chose 
to keep them both to estimate their effect on our results. As we show below, it is quite 
small. 

Finally, different diffusion schemes lead to different forms for the energy dependence 
of the diffusion coefficient and the reacceleration term. Several aspects of the diffusion 
process are treated in Schlickeiser 2002, and we considered three alternative possibili- 
ties : (i) Slab Alfvcn wave turbulence, with Ka(p) = Ko(37l s and K pp oc V^p 2 /Ka(j>), 
(ii) Isotropic fast magnetosonic wave turbulence, with Kp(p) = KqP 2 ~ s 7Z s and K^ p oc 
VaP 2 /3 1 - s Mv/V a )/K f (p) 1 and (iii) mixture of the two last cases, K M (p) = KqP^II 5 
and Kf p = Jf£. 

All results will be presented with the case (i), except in the specific discussion in 



Sec. £3j 

3.2. Summary: updates of Paper I's formulae 

The only changes with our previous study are 

— Eq. (19) of Paper I is replaced by 

,i , s / dE\ I dE\ /dE\ I dE\ 



Rcac 



where 



(1 + 2 ) 

(dE/dt) Rcac = [ --^K pp . (7) 



Eq. (A13) of Paper I (second order differential equation to solve) reads now - we use 
the same notations - 

A\ 'N* ; (0) = & 2h A ^ lo jE)N^0) - ^ Ajvi( )| (8) 
with 
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In our model, h roac = h, but it has to be kept in mind that a possible reinterprctation 
of V a is always possible (see Sec. |^) as long as h rea ,c <C L (this condition is necessary 
for the solution to be valid). 
— As regards the source spectra, two forms (hereafter type (a) and (b)) are used instead 
of Eq. (9) of Paper I 

a-Q(E) oc ~iT Q (10) 
b - Q(E) oc R- a (11) 

where R is the rigidity and a a universal slope of spectra for all nuclei heavier than 
helium. 



4. Runs and selection method 

The analysis presented here is the natural continuation of the work presented in Paper I. 
It is more general and it encompasses its results as a five-dimensional subset of the 
six-dimensional space scanned here. 

4-1- The six free parameters of the study 

The six parameters of this study are: the spectral index of sources a, the normalization 
Kq and spectral index 5 of the diffusion coefficient, the height of the diffusive halo L, 
the Galactic convective wind speed V c and the Alfvenic speed V a - They are included 
in our code as follows (sec Fig. [j] for a sketch of the procedure): for a given set of 
parameters, source abundances of all nuclei (i.e. primaries and mixed nuclei) are adjusted 
so that the propagated top of atmosphere fluxes agree with the data at 10.6 GeV/nuc 
(see Paper I). We remind that for B/C ratio, we checked that starting the evaluation 
of fluxes from Sulfur is sufficient (heavier nuclei do not contribute significantly to this 
ratio). Top of atmosphere fluxes are deduced from interstellar fluxes using the force field 
modulation scheme (sec Paper I and references therein). The resulting B/C spectrum 
is then compared to the data (see below) and a \ 2 1S computed for the chosen set of 
parameters. 

This procedure is very time consuming. Even when the location of \ 2 minima in the 
six-dimensional parameter space are known, more than 2.10 6 configurations are needed 
to have a good sampling of the regions of interest, for a given form of the source term 
energy dependence. 
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Fig. 1. Diagrammatic representation of the various steps of the propagation code 



4-2. x 2 criterion of goodness 



As in our previous analysis, we have computed the quantity 

2 _ ((B/@)i,exp ~ (B / model) 



\ 



E 



(12) 



where the sum runs over 26 experimental values from heao-3 (Engclmann et al., 1990) 
with energies ranging from 620 MeV/nuc to 35 GeV/nuc (as in Paper I). In general, if the 
experimental set-up is such that the measured (experimental) values differ from the "real" 
values by a quantity of zero mean (non biased) with a given probability distribution, 
then the value of x 2 gives a quantitative estimate of the probability that the model is 
appropriate to describe the data. However, this condition is probably not fulfilled for 
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HEAO-3, as for some measured quantity, the quoted errors (J^ exp are much smaller (e.g. 
oxygen fluxes) or much larger (e.g. sub-Fc/Fe ratio) than the dispersion of data itself. 
For this reason, it is meaningless to associate a likelihood to given \ 2 values. Instead, in 
Paper I we decided that models giving x 2 less than some value Xo were "good fits" while 
the others were "poor fits" . In this paper, no cut is applied and all the models, whatever 
the value of x 2 1 are shown in the figures. 

5. Results 

5.1. Subset 1: fixed measured spectral index 7 = a + <5 = 2.8 

In this section we present the results obtained for source spectra of the form Q(E) cx 
R~ a / (3 and diffusion coefficient K = Kof3TZ s . At sufficiently high energies, spallations 
and energetic changes are irrelevant and the measured fluxes can be considered as a mere 
result of acceleration and diffusion (see for example Maurin, Casse <fc Vangioni-Flam 
2002). In this case, the observed spectrum is proportional to 1Z~ 7 with 7 = a + S. In 
this section, we focus on the situation 7 = 2.8, corresponding to the spectral index of the 
measured Boron progenitor fluxes. Actually, Wiebel-Sooth, Biermann & Meyer (1998) 
analysed data from several experiments and derived smaller values. In Paper I, we found 
that the Oxygen flux measured by heao-3 would be more compatible with our diffusion 
model for a higher 7, namely 2.8 instead of 2.68. Anyway, we are more interested in the 
trends in the variation of other parameters for a fixed value of 7 than in this precise 
numerical value. The other cases will be treated in the following sections. 

We also set the halo thickness L to 6 kpc, leaving us with four free parameters (S, 
Kq, V a and V c ). All curves depicted in Fig. || correspond to one-dimensional cuts through 
the absolute x 2 minimum (for a given <5, the three different cuts justify the fact that we 
are located in a minimum). In the upper panel of Fig. ^, we plot the values of the x 2 as 
a function of Kq/L, for different values of S (and the corresponding a = 2.8 — <5). The 
best fits are obtained for S ~ 0.8-0.9, far from the Kolmogorov spectrum (5 = 1/3). We 
found a quite similar result in Paper I, where the same assumptions on 7 were made but 
with a different choice for the source spectrum, Q(E) = TZ~ a . The fit is best for values 
of the diffusion coefficient normalization Kq ~ 6 x 10~ 3 kpc 2 Myr~ x , yielding the value 
X 2 ~ 40 (giving x 2 ~ 1.8). For a Kolmogorov spectrum, the minimum x 2 is almost twice 
this value. Leaving aside any statistical interpretation of the analysis, we can observe 
that for greater 5, the minima of x 2 are obtained for smaller Kq/L (or Kq, L being set 
to 6 kpc) and versa- vice. This can be understood as at a sufficiently high energy Ethrcsh, 
diffusion is the sole remaining influential parameter and, for the flux to be unchanged 
with various 6, one need to satisfy roughly the relation Kq x -Efhresh ~ cte (this will also 
explain why type (a) and type (b) source spectra give similar Kq, see below). 
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Q(E) = R~ a / i 8 , 7 = ^ + 6 = 2.8 
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Fig. 2. Evolution of the x 2 value for various combination of parameters. All curves are 
for type (a) spectra (R~ a /f3) with 7 = 5 + a = 2.8 and the halo size is fixed to L = 6 kpc. 
Each curve shows one-dimensional cuts in the Kq/L (upper panel), V c (left lower panel) 
or V a (right lower panel) direction of the 3-dimcnsional \ 2 hyper-surface, fn the upper 
panel 5 is varied from 1.0 to 0.3 and, in the lower panels the same symbols to indicate 6 
are conserved. Each curve gives the absolute minimum for the parameter on the abscissa 
axis, L being fixed to 6 kpc (similar curves with slightly different minima are obtained 
for other L values). 

In the lower panels we present two cuts in the two other directions, namely in the 
V c and V a directions. The first one tells us that except for the special case 8 = 0.3 for 
which the \ 2 curve skips to null V c , B/C is fitted with V c between 10 and 20 km s _1 . 
The best \ 2 are for convective velocity around V c ~ 16-18 km s _1 . For V c < 15 km s _1 
(and 8 < 0.6) the goodness of the fit quickly decreases. We can see that when 8 is around 
0.4-0.3, the B/C ratio becomes very sensitive to the V c values. It appears that when 8 
is decreased, a good fit is maintained provided that V c is also lowered. This is possible 
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down to S ~ 0.4 for which the best value for V c is zero. For lower <5, the previous trade-off 
cannot be achieved (as V c must be positive for the galactic wind to be directed outwards) 
and no good fit is possible. 

The right panel shows the \ 2 curves as functions of the Alfven velocity. The mini- 
mization procedure always yields a V a far different from zero. Good fits are obtained for 
values of V a ~ 40-50 km s _1 . 

In each of the explored directions, the \ 2 curves are very narrow: the diffusion model 
leads to meaningful and intcrprctablc values for all the physical, free parameters. Similar 
results, with slightly different values for the minima, are obtained for the other values of 
L in the range 1 < L < 15 kpc. 




Fig. 3. Left panel: evolution of the best \ 2 value with Kq/L for various S (1.0 to 0.3, 
from left to right) at a fixed 7 = 2.8. Each curve correponds to a given halo size L from 
14 kpc to 2 kpc. Right panel: the same best \ 2 values are presented versus V c and V a . In 
both panels, empty circles correspond to type (a) spectra and stars to type (b) spectra. 



In Fig. g we present the results for the same analysis for different values of the halo 
thickness L and considering also the form (b) for the source spectra, i.e. Q(E) oc Rr a . 
The total spectral index 7 is still set to 2.8. The left panel reports the \ 2 as functions 
of Kq/L, for different values of S and L, and for both types of source spectra. We see 
that the choice (b) globally improves the fit, and the favoured range for S is now 6 > 0.4 



Maurin, Taillet & Donato: Source and diffusion spectral features. . . 11 

(whereas 5 > 0.7 for choice (a)). At fixed S and L, the absolute minima for both choices 
correspond to very similar values of Kq/L. We can also notice that type (a) spectra are, 
for the higher S, more sensitive to variations of L. 

In the right panels we show a cut in the V c -V a plane. For both type (a) and (b) 
spectra, (5=0.3 yields a null value for the convective wind. Type (a) spectra give a little 
bit higher V c . At fixed S, the variation of L has almost no effect on V c , while it is strongly 
correlated with the increase of V a . 

5.2. Subset 2: S = 0.6, new features from a variation 

In this section we discuss the results obtained when the index a is varied between 1.3 
and 2.5, S being set to a given value 6 = 0.6 which has been extensively used in the 
literature. 



Q(E) = R~%3 , 6 = 0.6 L=6 kpc 




0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4 0.42 




V c (km/s) V ft (km/s) 



Fig. 4. Same as in Fig. || (type (a) spectra, L = 6 kpc), but for a fixed 5 = 0.6. 
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Fig. H corresponds to the previous Fig. |[ In the left panel we observe that a large 
variation of the index a has a slight effect on the normalization of the diffusion coefficient 
Ko, which stays around an average value Kq/L ~ 0.0032 kpc Myr -1 for L = 6 kpc. 
Evolution of the absolute \ 2 minimum is also far less sensitive to a than S (see previous 
section). However, for a > 2.2 the fit to the data is poor and a global power 7 > 2.8 at 
S = 0.6 is excluded. 

The lower panels represent a cut in the V c and V a directions. We can observe that the 
minimization procedure always drives the minima towards convective velocities between 
12 and 16 km s _1 , the least V c being obtained for the smallest a. This range is again 
very narrow. Similarly, reacceleration is needed to fit data and the minima of the x 2 are 
obtained for V a between 55 and 75 km s" 1 . Towards this lower limit, \ 2 is high and the 
model cannot confidently reproduce observations. 

When 5 is fixed, we can conclude that a variation in the power of the type (a) source 
spectrum does not strongly act on the evolution of Kq, V c and also V a . This can be also 
easily understood : forgetting for a while energy gains and losses, we see from diffusion 
equation solutions (the same behavior occurs in leaky box models) that the source term 
can be factorized so that secondary to primary ratios finally do not depend on Q(E), i.e. 
are independent of a. Once again, the absolute minimum is identified by a steep \ 2 m 
these three directions. 
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Fig. 5. Same as in Fig. |^ but for a fixed 6 = 0.6. 
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In Fig. H we present the results for 8 = 0.6, and for both type (a) and (b) source 
spectra, to focus on the evolution of L and a. The left panel tells us that the evolution 
of the halo thickness from 2 to 14 kpc, at fixed a (in other words, at fixed 7 = a + 0.6) 
does not change the goodness of the fit. Only a slight modification in Kq/ L is required in 
order to recover the same B/C flux ratio. Type (b) source spectra reproduce quite well 
the data for all the explored parameter space. On the contrary, the better theoretically 
motivated type (a) spectra cannot reproduce observations for a > 2.2 if 8 = 0.6. Since 
at high energies the two source spectra are equivalent, we must conclude that it is the 
low energy part of B/C which is responsible for such a discrimination. 

The right panels show the absolute minima in the V c -Va plane. Both spectra require 
non null reaccclcration and convection. Even more so, the selected values reside in the 
narrow interval for V c , i.e. V c ~ 10 — 15 km s _1 and between 40 and 90 km s _1 for the 
Alfven velocity. 



5.3. Subset 3: a = 2.0, standard acceleration 

Fig. ^ describes the results of the analysis done assuming type (a) source spectra, with 
fixed index a = 2.0 and L = 6 kpc. A consensus seems to emerge in favor of values 
a sa 2.0 (see Drury et al., 2001 and references therein), close to the index given by 
primeval acceleration models, but any other value would be fine for the purpose of this 
section. In the upper panel 8 has been varied between 1.0 and 0.3, and the figure shows 
the evolution of the \ 2 with respect to Kq/ L. As in Fig. ^ and, at variance with Fig. [|, 
the minima correspond to Kq/L spanning over almost two orders of magnitude. It is the 
modification of the power-law in the diffusion coefficient - and not in the source spectrum 
- that significantly acts on Kq. Once again, the Kolmogorov spectrum is disfavoured: in 
this case it is obvious that the calculated flux ratio would be too hard. The best fits are 
obtained for 8 ~ 0.6 - 0.9. 

The lower panels show the cuts in the V c and V a directions. The left one tells us that 
for smaller 8, the preferred convective velocities are smaller (and the best \ 2 IS larger), 
down to 5 = 0.3 for which a no-convection model is prefered, with a bad \ 2 - The best 
fits are obtained for V c around 15-18 km s _1 . In the right panel we can notice, again, 
that only models with reaccclcration have been chosen by the minimization procedure. 
Lower 8 point to higher Kq/L and V a values and lower V c . The same trend is recovered 
in the other cases treated above. Reacceleration and convection act, in a certain sense, 
in competition, even if data always give preference to a combined effect rather than their 
absence. 

This trend (the smaller 8, the larger Kq, or equivalently Kq/L as L is constant in 
the above figures) was already mentioned in Sec. 5.1 , Actually, as we will see in Sec. ^, 
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Fig. 6. Same as in Fig. || (type (a) spectra, L — 6 kpc), but for a fixed a = 2.0. 

the correlation between Kq/L and V c is more properly explained by virtue of Eq. ( |l7| ) 
so that the evolution of V c is fixed by the evolution of the two other free parameters, 
i.e. Kq/L and S. As regards V a , it only appears in Eq. (||). A rough estimation can be 
inferred using power-laws K(E) cx E 5 and N ] {E) cx E~ ( - a+5 '> in Eqs. (@) and @: 

One finally obtains that the term for energetic redistributions evolves as (3 + <5)(1 + 
26)N(E) x Va/K(E) for a = 2. Hence, from the above argument, when S is decreased, 
Kq is adjusted so that K{E) and N(E) remain grossly the same. However, for the above 
expression to be constant, V a must be increased; this is the trend we observe. 

In Fig. ^ we show the effect of varying the halo thickness when the source spectral 
index is fixed to 2.0 and all the other free parameters are scanned. Again, type (b) spectra 
reproduce better the data. When L is varied between 14 and 2 kpc, this may modify the 
chosen Kq/L by a factor of two. The right panels tell us that the influence of L on V a 
is to double its value when L is varied from its minimum to its maximum value. On 
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Best x 2 for various 8 and L combinations (fixed a = 2.0) 
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Fig. 7. Same as in Fig. || but for a fixed a = 2.0. 

the contrary, the effect on V c is almost null. The situation for V c and V a is very similar 
to the one discussed in the two above cases, when 7 and then a were fixed. Indeed, 
looking carefully at the above figures, we recover the same effect also for K /L, at fixed 
a + 5. Again, the behaviour of V c can be understood but cannot be simply explained. 
Conversely, neglecting V c in the asymptotical formula , one can see that when L increases, 
Kq must increase (as can be checked in the left panel). Moreover, it can be seen from the 
form of K pp that V a increases as the square root of Kq when S is fixed (see right lower 
panel). 

5.4- The whole set: final results 

We know present the result of the full analysis, in which all the parameters are varied. 
Fig. U shows the evolution of the x 2 m the 5 and 7 = S + a plane for different values of 
L. We can see that, at fixed type (a) or (b) spectra, a change in the halo height L has 
almost no effect on the best \ 2 surface. Generally, high values for 5 are preferred and, 
a Kolmogorov regime for the spatial diffusion coefficient is strongly disfavoured over all 
the parameter space. More precisely, type (b) spectra point towards a band defined by 
S ~ 0.8 in the £-7 plane, whereas the type (a) spectra gives the additional constraint 
7 < 2.8 (sec Fig. |H). 
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Fig. 8. Best x 2 values for various L (2, 6 and 10 kpc) in the plane 8 — 7. Left histograms 
are type (a) spectra and right histograms type (b). Notice that for right histograms, only 
the upper figure displays the values <5 = 0.3 and S = 0.35. They have been omitted in the 
two remaining figures to gain contrast (for any L, these configurations have x 2 ^ 100). 
Assuming L = 6 kpc, type (a) source spectra give a best value Xbest = 1^.8 for a = 1.65 
and S — 0.85 whereas type (b) gives Xbost = 14.6 for a — 1.95 and 6 = 0.85. These were 
obtained with 26 data points. 
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In Fig. |^ we show the preferred values of the three remaining diffusion parameters 
Ko, V c and V a , for each best x 2 m the S-j plane, when L has been fixed to 6 kpc. The 
two upper panels show that the evolution of a does not affect Kg. On the other hand, 
as already noticed, we clearly see the (anti)correlation between the two parameters Kq 
and S entering the diffusion coefficient formula, giving the same normalization at high 
energy (K x Ef hrcsh sa cbe). Almost the same numbers are obtained for type (a) and (b) 
spectra. Kq spans between 0.003 and 0.1 kpc 2 Myr -1 . We will discuss in the following 
sections how these results can be compared to the literature. 

The middle panels show the values for the convective velocity. Only very few con- 
figurations include V c = 0, always when 5 = 0.3, for both types of source spectra. The 
value of V c increases with i5. For type (a) spectra, increasing 7 and S at the same time 
makes V c change its trend. As remarked previously, the effect of Galactic wind is more 
subtle since it acts at intermediate energies and is correlated with all the other diffusion 
parameters through the numerous terms of the diffusion equation. 

The lowest two panels show the influence of V a . We recover a correlation similar to 
the one discussed for Kq (see Eq. |l3|). The Alfven velocity doubles from 5 = 1.0 to 0.3, 
whereas it is almost unchanged by a variation in the parameter 7 (or equivalently a). 

All the three analysed parameters (i.e. Kq, V c and V a ) behave very similarly with 
respect to a change in the source spectrum from type (a) to type (b). It can be explained 



as the influence on the primary and secondary fluxes can be factored out (see Sec. 5.2) if 
energy changes are discarded (their effect is actually small on the derived parameters). 
Existing data on B/C do not allow us to discriminate clearly between these two shapes 
for the acceleration spectrum. This goal could be reached by means of better data not 
only for B/C but also for primary nuclei (Donato et al., in preparation). 

5.5. Other diffusion schemes 



As discussed in Sec. 3.1, we tested three different diffusion schemes, with three different 
forms for the diffusion coefficient. Most results are basically insensitive to the choice of 
this form. In particular, the figures corresponding to Fig. |9| are almost identical to the 
case presented above, so that they will not be reproduced here. Figure |l^ displays the 
X 2 as a function of S and 7. The values of x 2 are slightly different in the three cases, but 
the general trend is the same, and all the previous conclusions still apply. 

5.6. Sub-Fe/Fe ratio 

In an ideal situation in which we had very good and consistent data on B/C and sub- 
Fe/Fe ratios, the best attitude would be to make a statistical analysis of the combined 
set of data. Unfortunately, this is not currently the case. We consider two ways to extract 
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Fig. 9. From top to bottom: for each 
corresponding values of log(-Ko), V c and 




6 



jst x 2 m the plane J — 7 (L = 6 kpc), the 
j are plotted for both source spectrum types. 



information from the Sub-Fc/Fe data. First, as a check, we compare the sub-Fe/Fe ratio 
predicted by our model - using the parameters derived from our above B/C analysis - 
with data from the same experiment. Second, we search directly the minimum %p e of the 



sub-Fe/Fe ratio, with no prior coming from B/C. As previously emphasized (see Sec. 4.2), 
this procedure is more hazardous since the statistical significance of the sub-Fe/Fe data 
is far from clear. 



5.6.1. Using B/C-induced parameters to derive Xb/c^Fc 

For each set of diffusion parameters giving a good fit to the observed B/C ratio, the 
sub-Fe/Fe ratio can be computed and compared to the values measured by heao-3. This 
is not as straightforward as in the B/C case because although Sc, Ti and V - that enter in 
the sub-Fe group (as combined in data here) - are pure secondaries, some of the species 
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K (S 2 ' 6 R s 




K (S'- s R s 

Fig. 10. Best x 2 values, in the plane 5 — -f, for the three different forms of the diffu- 
sion coefficient and reacceleration terms (i) Slab Alfven wave turbulence, with Ka{p) — 
Ko/37Z s and K pp oc V\p 2 / K a{p) , (ii) Isotropic fast magnetosonic wave turbulence, with 
K F [p) = K P 2 - 5 K S and K^ p oc Vffifi 1 - 5 \n(v/V A )/K F (p), and (iii) mixture of the two 
last cases, K M (p) = KqP 1 '^ 5 and = KfL. 



20 Maurin, Taillet & Donato: Source and diffusion spectral features. . . 

intermediate between sub-Fc and Fe. contributing to the sub-Fe flux, are mixed species 
(i.e. Cr, Mn). As a consequence, all the primary contributions were adjusted so as to 
reproduce the sub-Fe/Fe ratio at 3.35 GeV/nuc. The sub-Fe/Fe spectra are not steep 
enough at high energy, so that normalization at 10.6 GeV (i.e. as for B/C) would have 
led to less good fits. We emphasize that to perform this normalization of secondary-to- 
primary is equivalent to making an assumption about the elemental composition of the 
sources, which is usually deduced from secondary-to-primary ratios. A different choice 
would slightly shift the normalization of sub-Fe/Fe ratio without affecting much our 
conclusions. 

Fig. [ll] displays the Xb/c^Fo vamGS obtained when the diffusion parameters giving a 
good fit to B/C are used to compute the sub-Fe/Fe ratio , for each value of a and 5 (for 
type (a) spectra and L — 6 kpc, although the results for type (b) and/or different L are 
quite similar). This surface is very similar to the surface obtained with B/C, pointing 
towards high values of S (compare to Fig. |J). 

5.6.2. Looking for XF e 

We now consider a full sub-Fe/Fe analysis (i.e. the parameters minimizing Xp e are looked 
for) but we emphasize that the results given here are from our point of view far less robust 
than those obtained from B/C. As a consequence, conclusions of this section have to be 
taken only as possible trends. Several points can be underlined from Fig. [l^: (i) as for the 
B/C case, the best Xfc ^ s obtained for type (b) spectra, (ii) the general behavior of Ko, 
V c and to a less extent V a is mostly the same as for B/C. (iii) the type (b) spectra yield 
propagation parameters which are closer to B/C's, as we can see from V c values; (iv) 
finally, consistency with B/C analysis would be better obtained for S pointing towards 
0.6-0.7. 

5. 7. Additional insight from visual comparison of our model to data 

Typical spectra (modulated at $ = 500 MV) are shown in Fig. [l^, for different values of 
the parameters a and 5, along with the data points from HEAO-3 (Engclmann et al., 1990) 
and balloon flights (Dwyer & Meyer, 1987). Three low-energy data points, from HET on 
Ulysses (Duvernois & Thayer, 1996), HKH on ISEE-3 (Leske, 1993) and Voyager (Webber, 
Lukasiak & McDonald, 2002) are also shown; they all have about the same modulation 
parameter, i.e. $ « 500 MV. The ACE points ($ « 750 MV) are also displayed (Davis et 
al., 2002). 

All the models displayed give similar spectra, which would be difficult to sort by eye. 
This may explain why some of these models (e.g. those with S = 0.3) are retained in 
other studies. The main features are (i) the influence of S on the high energy behaviour - 
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Fig. 11. Values of Xb/c^Fc obtained by applying the diffusion parameters - type (a) 
source spectra and L = 6 kpc - giving the best fit to B/C, for each a and 6, to sub- 
Fe/Fe. The X B / c ^ Fc are computed with HEAO-3 data points. 

a good discrimination between these models would be provided by precise measurements 
around 100 GeV/nuc - and (ii) the type (a) source spectra are steeper than type (b) at 
low energy. 



6. Comparison with other works 

Some of our configurations can be compared to those previously found in similar models. 
In particular, to compare the Alfven speed from one paper to another, we have to be 
sure that all V a used denote the same quantity. 

To compare the reacceleration terms employed, we retain only the spallation term 
and the highest order derivative in energy in the diffusion equation, giving 
2hn lSM va 3 N J (0) = 2hf3 2 K pp d f^ 0) . (14) 
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Q(E) = R"//S (L . 6kpc) Q(E) = R 




Fig. 12. From top to bottom: best Xfc an< ^ f° r eacn best x| e in the S — 7 plane (L = 6 
kpc), the corresponding values of log(-Ko), V c and 14 are plotted for both source spectrum 
types. 



We have supposed that both phenomena occur only in the thin disc h <C L and, in the 
above equation, the reacceleration zone height equals the spallative zone height. If it is 
not the case, we have to correct the previous relation by a multiplying factor h iea , c /h. 
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Fig. 13. The B/C and sub-Fe/Fe spectra (modulated at $ = 500 MV) for several sets 
of parameters (giving the best fit to B/C for these values) are displayed, along with 
experimental data from heao-3 (Engclmann ct al., f 990), ballon flights (Dwyer & Meyer, 
1987), HET on Ulysses (Duvernois & Thayer, 1996), hkh on ISEE-3 (Leske, 1993) and 
Voyager (Webber, Lukasiak & McDonald, 2002). Note that ACE data (Davis et al., 2002) 
correspond to a modulation parameter $ « 750 MV. 

Actually, V a is "fixed" through the choice of K pp K(E). As underlined in Sec. |3.l| , this 
paper now follows the requisites of minimal reacceleration models (see Tab. [j], last line). 

Once this h reac /h g!ijS rescaling - that differs from one paper to another - is taken into 
account, a comparison is possible between models if a minimal resemblance exists between 
the other input parameters, i.e. same 5, a (plus same form of the source spectrum) and 
halo size L; Tab. [I] shows the value adopted for these parameters in two recent studies. 

6.1. Maurin et al. (2001) - Paper I 

The results are expected to be slightly different from our previous study as the compo- 
nents have been modified. First, V a has a different interpretation in the two studies (see 
Tab. |l|, first column). As underlined above - remembering that in Paper 1 the diffusion 
coefficients scaled as K PP K(E) = (2/9) x p 2 V£ -, the Alfven speed value from Paper 1 
(V r Q Papor : ) has to be rescaled into V^ >aper standard (i.e. as the standard convention used 
in this work and others) through the relation 

/ 6 

T/Paper I r/Papcr I, Standard /l K\ 

Va ~ Va V S(4-5*)(4-6) ■ (15) 

Second, the equation describing diffusion in energy has been modified and, the values 

of K , V c and V a that give the best fit to B/C data for a given 8 must change at some 

level. Notice that in Paper I we used a source term corresponding to type (b) spectra 
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Table 1. Main characteristics of various diffusion models 



Maurin et al. / This work 
(2001) / (2002) 



Seo & Ptuskin / Jones et al. 
(1994) / (2001) 



Moskalenko et al. 
(2002) 



Thin disc h (pc) 
Halo size L (kpc) 
2/f 

Surface mass density* 
K PP K{E) 



h = 100 pc 



2h le 



■lit 



10" 3 g cm" 2 



h = 200 pc 
L — 3 kpc 

2/f 

reac — 2L/3 f 
2.0 x 10" 3 / 2.4 x 10" 3 - —~ 2 



Gas distribution 
L = 4 kpc 

2/ircac — 27j 



g cm 



1.6 x 10~ 3 



g cm 



2 p v i 



I 



4 P V a 



4 P V a 

3,5(4-5' A )(4-6)w 



4 P V a 

M(4-^)(4-<5)iij 



* Jones et al. use the same set of equation and parameters as Seo and Ptuskin for the stochastic 
reacceleration model. Consequently, it seems that h a defined in Jones et al. is the half-height of 
the reacceleration zone, contrarily to what is depicted in their Fig. 4. 

* The surface mass density is defined as /i = 2hp where p is the matter density in the thin disk. 

(see Eq. |ll|), with 7 of each species that were set to their measured value (see details in 
Paper I); this corresponds roughly to 7 w 2.7 for all boron progenitors. 

However, we find that the conclusions raised in Paper I, in particular the behaviors 
reflected in Figs. 7 and 8 of Paper I, are basically unchanged (it is not straightforward to 
compare with present figures, but the careful reader can check this result using the above 
scaling relation and the corresponding parameter combinations). To be more precise, it 
appears that Kq/ L does not significantly change (for example, for 6 = 0.6 and L = 2 kpc, 
we still have K$/L ~ 0.004 kpc Myr -1 , see Fig. || left panel - this paper - and Fig. 7 
of Paper I). As regards the galactic convective wind, V c is shifted towards higher values, 
whereas the V a /^/Ko range remains roughly unchanged. 

This can be easily understood: the additional term - comparable to a first order gain 
in energy, see Eq. (^) - has to be balanced to keep the fit good. This balance is ensured by 
enhanced adiabatic losses, i.e. bigger V c . Other parameters are only very slightly affected 
by this new balance. 



6.2. Jones et al. (2001), Moskalenko et al. (2002). 

Moskalenko et al. (2002) (hereafter Mos02) use a description more refined than ours 
because they include a realistic gas distribution. Jones et al. (2001) (hereafter JonOl) 
take advantage of an equivalent description in terms of a leaky box formalism (use of 
a phcnomenological diffusion coefficient) - for both wind model (no reacceleration) and 
minimal reacceleration model (no wind) - to solve the diffusion equation. 

Let us make a few comments at the qualitative level. First, starting with Mos02's 
models, we can note that convection has always been disfavored by these authors. For 
example, in their first paper of a series (Strong & Moskalenko, 1998), a gradient of 
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convection greater than 7 km s _1 kpc -1 was excluded. We notice that this result was 
not very convincing since it is clear from an examination of their figures that none of the 
models they proposed gave good fits to B/C data. Thanks to many updates in their code, 
their fits were greatly improved (Strong & Moskalcnko, 2001; Moskalenko ct al., 2002) 
but if it is now qualitatively good, it is hard to say how good it is since no quantitative 
criterion is furnished. Anyway, our Fig. |^ allows us to understand why convection is 
disfavored in such models. Actually, if 6 ~ 0.3 - as in the Kolmogorov diffusion slope 
hypothesis 6 = 1/3 -, we see that for such a configuration, the best fits are obtained for 
14 ~ km s -1 . 

Similar comments apply to JonOl's models. Given a Kolmogorov spectral index for 
the diffusion coefficient, their combined fit to B/C plus sub-Fe/Fe data is not entirely 
satisfactory. It improves for higher values of 8 and in the convective model (they do not 
include reaccclcration in this model), their best fit being obtained for 8 = 0.74. As the 
authors emphasized, the search in parameter space was not automated and they cannot 
guarantee that their best fit is the absolute best fit. Actually, the sub-Fe/Fe contribution 
to the x 2 value has to be taken with care. First, the error bars are not estimated well 



enough to give a statistical meaning for % 2 values (see Sec. 4.2) and a different weight 
should be considered for B/C and sub-Fe/Fe. Second, if the best parameters extracted 
from B/C data reproduce formally the same \ 2 surface when applied to the evaluation 
of sub-Fe/Fe (see Fig. [n]), the direct search for the parameters minimizing \ 2 f° r the 
same sub-Fe/Fe data gives constraints that are much weaker (see Fig. [l2]). Thus, any 
conclusion including this ratio is from our point of view far less robust. 



6.3. Quantitative comparison, interpretation of Kg, V c and V a 
6.3.1. Justification of the differences between models 

Tabs. H and || give the results of Mos02 and JonOl - without any rescaling of any param- 
eters - compared to what is obtained here; only a few models are displayed. 

Tab. H shows Kq, V c and V a for 8 = 0.30 and a = 2.40. Taking the first three lines at 
face value, our values of Kq and V a are very different from the others and our model seems 
to have a problem. However, the matter disk properties (height and surface density) are 
different in these models. To be able to compare, we set these quantities to the values 
given in JonOl and the resulting parameters are shown in the last line of Tab. |^. 

Actually, we know that in diffusion models, the behavior is driven by the location of 
the closer edge, leading to a preferred escape on this side. With L = 3 kpc, our three- 
dimensional model should behave as the two-dimensional model with infinite extension 
in the r direction of JonOl. This hypothesis can be validated if one takes their Eq. (3.6). 
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Table 2. Diffusion parameters obtained in models with S = 0.30, a = 2.40 for pure 
power-law source spectra. 



L 


h 


H x 10~ 3 


^reac 


K 


v c 


V a 


xl 




Ref. 


(kpc) 


(kpc) 


(g cm -2 ) 


(kpc) 


(kpc 2 Myr^ 1 ) 


(km s^ 1 ) 


(km s" 1 ) 








4. 


n(r) 


1.6 


4.0 


~ 0.201 


0. 


30. 


Good 




Mos02 § 


3. 


0.2 


2.4 


1.0 


~ 0.196 


0. 


40. 


1.8 




JonOl 1 " 


3. 


0.1 


1.0 


0.1 


~ 0.0535 


0. 


105.8 


4.2 


(Figs 


| and ||, this paper)* 


3. 


0.2 


2.4 


1.0 


~ 0.127 


0. 


47.3 


4.4 


This work* 



5 For this model, the exact values are 8 = 0.33, a = 2.43. 

* This model give best fit to flux using a slightly modified form for the source; Q oc 
{R/2)- 2 ] 1 ' 2 . 

* Corresponds to the best fit for the presented L, a and 8 value. 



For the pure diffusion model (reacceleration and convection are discarded), one has a 
simple relation between [i, L and Kq through an equivalent leaky box grammage 

Xdtf = 2^ (16) 
A direct application of this result to our model with h = 100 pc (third line of Tab. ||) 

using the scaling fj, — > 2.4 x /i, leads to a rescaling Kq — > 2.4 x Kq, consistent with results 

of the fourth line. 

A similar expression may be obtained in the presence of galactic wind : in the wind 



model (their Eq. 4.6), one has 



X 



\x v 
2V C 



1 — cxp 



-V C L 



K„R S 

Applied to the second line of Tab. H, this gives V c 



(17) 



2.4 x V c , leading in turn to Ko 
2.4 x Kq 7 also in very good agreement with the direct output of our code. 

Table 3. Diffusion parameters obtained in models 5 = 0.74, a = 2.35. 



L 


h 


H x 10~ :i 


^reac 


K 


V c 


V a 


xl 


Ref. 


(kpc) 


(kpc) 


(g cm" 2 ) 


(kpc) 


(kpc 2 Myr" 1 ) 


(km/s) 


(km/s) 






3. 


0.2 


2.4 


1.0 


~ 0.024 


29. 


0. 


1.5 


JonOl 


3. 


0.1 


1.0 


0.1 


~ 0.0056 


15.5 


35.3 


3.0 


This work* 


3. 


0.2 


2.4 


0.1 


~ 0.0134 


36.5 


26.5 


3.1 





Corresponds to the best fit for the presented L, a and 8 value. 



Even with this fi rescaling, the diffusion coefficients obtained by the different authors 
quoted above are still not fully compatible. Another possible effect, namely the spatial 
distribution of cosmic ray sources, is now investigated. We note that in our model, the 
radial distribution of sources q(r) follows the distribution of supernovae and pulsar rem- 
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nants. The choice of this distribution has an effect on B/C spectra and on the parameters 
giving the best fits. If we use a constant source distribution q(r) = cte with V a set to to 
follow JonOl, we find that Kq is enhanced by about 10%. We checked that it is also the 
case for results presented in Tab. |[ Hence, it appears that results for 8 = 0.74 of JonOl, 
though slightly different, arc not in conflict with ours. As regards 8 = 0.3 and Mos2, 
using the scaling relation ( |l6| ) along with a 10% decrease of Kq for JonOl, we obtain 
respectively Kq = 0.226 (Mos02), 0.176 (JonOl), 0.127 (this paper) kpc 2 Myr" 1 . Thus 
there is some difference between Mos02 and JonOl, which is not obvious when the values 
taken naively from Tab. || are compared. These discrepencies could have several origins: 
treatment of cross-sections (we checked that total and spallative cross sections - taking 
into account ghost nuclei, see Paper I - are compatible with recent data, e.g. Korejwo et 
al., 2002), average surface density in Mos02 that is probably not exactly 1.6, choice of 
data and fits for JonOl that differ from ours (some of the point they used are significantly 
lower than heao-3's). Finally, the fact that we scan the whole parameter space can make 
a difference from manual search. To conclude, results are qualitatively similar, but a few 
quantitative differences remain. The intrinsic complications and subtleties of the various 
propagation codes make it difficult to go further in the analysis of these differences. 

6.3.2. Meaning of K 

The normalization Kq gives a measure of the efficiency of the diffusion process at a 
given energy. Its value can be predicted if (i) a good modelling of charged particles in a 
stochastic magnetic field and (ii) a good description of the actual spatial structure of this 
magnetic field, were available. It is not the case and the precise value of Kq is of little 
interest. Moreover, the presence of effects other than pure diffusion can be mimicked, at 
least to some extent, by a change in Kq. Eq. ( |l7| ) gives a whole class of parameters giving 
the same results and can be used to extract an effective value of Kq taking into account 
the effect of the size of the halo L and wind V c . This also explains the great range of 
values that can be found in the literature. 

This relation shows that there is also an indeterminacy of the absolute density of the 
model, because as long as h x nigM is constant, the grammage Xdif is also constant. 
Fortunately, a realistic distribution of gas can be deduced by more direct observational 
methods, so that a definite value of uism can be used. 

6.3.3. Galactic convective wind V c 

We note that in our model, Galactic wind is perpendicular to the disc plane and is 
constant with z. Actually, the exact form of galactic winds is not known. From a 
self-consistent analytical description including magnctohydrodynamic calculations of the 
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galactic wind flow, cosmic-ray pressure and the thermal gas in a rotating galaxy, Ptuskin 
et al. (1997) (see also references therein) find a wind increasing linearly with z up to 
z ~ f5 kpc, with a z = value of about 22.5 km s -1 . Following a completely different 
approach, Soutoul & Ptuskin (2001) extract the velocity form able to reproduce data from 
a one-dimensional diffusion/convection model. They obtain a decrease from 35 km s" 1 to 
12 km s _1 for z ranging from 40 pc to 1 kpc followed by an increase to 20 km s _1 at about 
3 kpc. For reference, our values for the best fits correspond to about 15 km s _1 . The 
difficulty to compare constant wind values to other z-dependences is related to the fact 
that cosmic rays do not spend the same amount of time at all z, so that there cannot be a 
simple correspondence (see also next section) from one model to another. As a result, all 
the above-mentioned models are formally different, with different inputs (spectral index, 
diffusion slope). Nevertheless, their values are roughly compatible, Ptuskin et al's model 
providing the grounds for a physical motivation for this wind. However, an even more 
complicated form of the Galactic wind could be relevant for a global description of the 
Galaxy (see Breitschwerdt et al., 2002). 



6.3.4. Interpretation of the Alfvenic speed V a 

Above, we gave some elements to compare V a values from various works. Actually, sec- 
ondary to primary ratios are not determined directly by the Alfven speed in the inter- 
stellar medium, but rather by an effective value: 
fh V cS 

"■-^¥4 (18) 

First, the parameter ui characterizes the level of turbulence and is often set to 1 (Seo & 
Ptuskin, 1994). Our model, as others, uses 

!1 if Z < /ireaci 
(19) 
otherwise; 

as a crude approximation of the more complex reality. 

Second, the total rate of reacceleration (at least in a first approximation, see discussion 
below) is given by a convolution of the time spend in the reacceleration zone and the 
corresponding true Alfven speed in this zone. There is a direct analogy with the case of 
spallations and the determination of the true density in the disc, as discussed above. The 
problem is still somewhat different, as there are no direct observational clues about the 
size of the reacceleration zone, or said differently, about u>(z). This leads to a degeneracy 
in h/h 1CBuC that holds as far as h rcac -C L, due to the structure of the equations in the 
thin disc approximation. For example, a model such as Strong et al's that uses h rcac = L 
cannot be simply scaled to ours. A cosmic ray undergoing reacceleration at a certain 
height z has a finite probability of escaping before it reaches Earth, this probability 
being greater for greater z. As a result, the total reacceleration undergone by a cosmic 
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ray is actually not a simple convolution of the reacceleration zone times the Alfven speed 
in this zone, but rather should be an average along z taking into account the above- 
mentioned probability (in principle, this remark also holds for the gaseous disk, though 
the latter is known to be very thin, ~ a few hundreds of pc). 

To conclude, there are basically three steps associated with three levels of approxima- 
tions to go from the V a deduced from cosmic ray analysis to the physical quantity. First, 
if u)(z) is approximatively constant with z, how large is the reacceleration zone height? 
The second level is related to the possibility that u>(z) strongly depends on z in a large 
reacceleration zone. If it is too large, the link with the phcnomenologically equivalent 
quantity in a thin zone is related to the vertical occupation of cosmic rays. However, 
this latter possibility seems to be unfavoured by MHD simulations (see Ptuskin et al., 
1997). Finally, with the above parallel between interpretation of /i and V a , we see how 
misleading it is to obtain precise physical quantities from our simple model, since there 
is no one-to-one correspondence between reality and simplified models. This discussion 
shows that even if the actual derived Alfven speeds are consistent with what is expected 
from "direct" observation (~ 10 — 30 km s _1 ), the cosmic ray studies would certainly 
be not very helpful in providing physical quantities better than a factor of two. If we 
reverse the reasoning and retain our best models with L = 6 kpc, we could conclude that 
y/h/hreac must be ~ 4 in order to give realistic values for V a (with evident a priori about 

6.3.5. The evolution of propagation models 

As suggested by the previous discussion, there are several propagation schemes, each 
associated with numerous configurations, that are able to explain the B/C data. Thus, 
the discussion should not be about the correctness of all these models (leaky boxes, two 
or three-dimensional diffusion models and their inner degeneracy) , but rather about their 
domain of validity. As a matter of fact, they are all equivalent, as far as stable cosmic 
rays around GeV/nuc energies are considered. 

Starting with the leaky box; it has been shown more than thirty years ago (Jones, 
1970) that the concept of "leakage-lifetime" was appropriate for the charged nuclei con- 
sidered here (see also Jones et al., 1989), even if it broke down for e~ (all orders in 
the development in "leakage" eigenmodes contribute because of synchrotron or inverse 
Compton losses) and for radioactive nuclei (Prishchep & Ptuskin, 1975). The leaky box, 
due to its simplicity, is very well suited for the extraction of source abundances (elemen- 
tal as well as isotopic). It can also be used for secondary antiproton production, since 
the same processes as for secondary stable nuclei are at work. However, as emphasized 
in Paper II, leaky box models are not able to predict any primary contribution in the 
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antiproton signal, since it requires the knowledge of the spatial distribution of primary 
progenitors. Considering a possible extension of leaky box models for stable charged nu- 
clei to high energy (~ PeV), it has been demonstrated in Maurin et al. (2002) that they 
are to a good approximation sufficient to describe the evolution of cosmic rays. Last, it 
is well known that leaky box parameters are just phenomenological with only a distant 
connection to physical quantities. 

This was further realized by Jones (1978, 1979) who first remarked that the phe- 
nomenological behavior of the escape length at low energy could be due to the presence 
of a Galactic wind. Jones et al. (2001) investigated further this idea and generated several 
equivalent phenomenological escape lengths from several possible physical configurations 
of a one-dimensional diffusion model. The relation between one-dimensional models and 
leaky box models is thus firmly established and very well understood. Moreover, this 
relation elucidates some of the physical content of leaky box models. Now if one wishes 
to overcome the inherent limitations of these models and say, to compute some primary 
antiproton component, one has to go through a three-dimensional model. It is likely that 
these models can also be related to the Jones et al. models (see Taillet & Maurin, in 
preparation). Several arguments used in the previous sections illustrate this view, but 
this occurs at least if the halo size is small compared to radial extension of the Galaxy. 

In the semi-analytical two-zone model used here, it is possible to evaluate the primary 
antiproton component (see Barrau et al., 2001) and to take into account radioactive 
species, even in the presence of a local very underdense bubble (see Donato et al., 2002 
for details). Our model fails to consider species such as e~ and e + , since the latter suffer 
from large energetic losses in the halo so that no simple semi-analytical approach can be 
used. The parameters extracted from these models are much easier to interpret in terms 
of physical quantities. 

Most of the limitations mentionned above are overcome by Strong et al's models. 
In this fully numerical model, all cosmic ray species can be computed self-consistcntly 
with the same propagation parameters. The main difference with our model is that a 
more realistic matter distribution is used instead of a thin homogeneous disk. They also 
consider that reacceleration occurs in the whole diffusion halo, which in our opinion is an 
approximation no more justified than the fact to confine it in a thin disk (see discussion 



in Sec. 6.3.4). Considering the gas distribution, both models are equally predictive for 
the charged nuclei (including antiprotons, see Fig. 9 of Paper II). On the one hand, our 
approach is better suited to scan the whole parameter space as we did in Paper I and in 
this paper. On the other hand, the Strong et al. models can check the consistency of e~, 
e + and 7 with observations, and can include whatever deviation from ideal cases for Kq, 
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V a , V c and more generally for any ingredient that enters in the description of propagation 
models. 

To conclude about models and their use, Jones et al.'s approach is probably the best 
and simpliest way to understand how physical parameters affect the propagated flux. Our 
model is very well suited for a consistent evaluation of all charged nuclei and extraction of 
propagation parameters; furthermore it is an intermediate step where general behaviors 
can still be analytically explored (Taillet & Maurin, in preparation). In the Strong et 
al. model, all fine effects can be studied and modelled, with the counterpart that the 
numerical approach makes the physical intuition of the results less straightforward. In its 
present form, Strong et al's model can be viewed as a fully numerical version of ours, so 
that their behaviors are very close. This discussion could leave the reader with a feeling 
that apart from these different modellings left to personal taste, galactic propagation 
phenomena are well understood. It is surely not the case! Even if all these models are 
equivalent to describe the local observations of charged cosmic rays, they lead to very 
different conclusions and interpretations when the spatial variation of the cosmic ray 
density is considered. As an illustration of the poor current understanding of this global 
aspect, we mention the ever-lasting problem of the gamma ray excess about 1 GeV 
towards the Galactic center or the too flat radial 7-ray distribution observed in the disc 
(see Breitschwerdt et al., 2002). 

7. Conclusion 

Forgetting for a while some of our theoretical a priori about the diffusion power spectrum, 
a new picture of cosmic ray propagation seems to emerge, motivated by the B/C analysis. 
In this new picture, high values for the diffusion coefficient spectral index (S > 0.6 — 0.7) 
and source spectral indices a ~ 2.0 are favored. This latter result is rather satisfactory: 
as emphasized in a recent working group report on SNR shocks (Drury et al., 2001), 
even if nonlinear acceleration models do not produce precise power-law spectra [...] the 
effective differential energy spectral index is close to 2.0. Furthermore, as pointed out in 
a series of papers by Vainio and Schlickeiser (2001, and see references therein), diffusive 
shock wave acceleration naturally yields smaller values of a if the correct scattering cen- 
ter compression ratio is used instead of the gas compression ratio. " 
This trend should be carefully analysed and discussed in the light of measured differential 
fluxes, in order to confirm or point out the possible inconsistencies in the current prop- 
agation treatments (see companion paper, Donato, Maurin & Taillet, in preparation). 
Briefly, the major arguments against large S come from anisotropy measurements at high 
energy and from theoretical preference for Kolmogorov-like turbulence spectra. However, 
Ptuskin et al. (1997) in their self-consistent analytical propagation model including gas, 
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cosmic ray and magnetic field - derived 5 ~ 0.55, a ~ 2.1 and argue that the observed 
anisotropy could be as well due to a particularity of the local structure of the Galactic 
magnetic field. Theoretical objections against too high values of S are probably more 
robust. 

For the rest, the conclusions of this paper can be summarized as follows : (i) we 
performed for the first time a full analysis of diffusion/convection/reacceleration models 
in the whole 6- dimensional parameter space (a, <5, Ko, L, V a , V c ), and S ~ 0.7 — 0.9; the 
values a ~ 2.0 are preferred; (ii) this preference holds whatever the specific form of the 
spectrum at low energy; the numerical values of the other parameters are also only slightly 
modified by this low energy dependence even though deviation from a power-law at low 
energy is preferred. The study of fluxes should give a more definite answer; (iii) Kq scales 
logarithmically with S and models with small halos tend to one-dimensional models with 
a simple relation between //, Kq, L and V c (see also Taillet & Maurin, in preparation); 
(iv) several existing models are compared and the qualitative and quantitative differences 
between them are studied and partially explained. 
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